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ABSTRACT 



Supersonic flow past oscillating flat plate cascades 
with supersonic leading-edge locus is analysed using a lin- 
earized method of characteristics valid for arbitrary fre- 
quencies and an elementary analytical theory valid only for 
low frequencies of oscillation. These two methods are ex- 
tensions of previous work by Teipel and Sauer for the single 
airfoil in an unbounded supersonic flow to the case of air- 
foils oscillating in cascade. Included is the determination 
of pressure distribution and both a two-degr ee-of -freedom 
(bending and torsion) flutter analysis and a single-degr ee- 
of-freedom (torsion) flutter analysis. Numerically deter- 
mined flutter boundaries are presented for various primary 
parameters such as, Mach number, solidity, stagger angle, 
density ratio, structural damping coefficient, and elastic 
axis position. Also, results are presented for the related 
problem of supersonic wind tunnel interference (including 
the effect of tunnel porosity) and airfoil-airfoil interfer- 
ence . 
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I. INTRODUCTION 



In order to design an effective supersonic compressor, an 
analysis of the complex unsteady flow phenomena in this speed 
regime is necessary. Since the infinite cascade has long been 
used by the engineer to model two-dimensional compressor 
flows, the study of oscillating supersonic cascades is import- 
ant in predicting the flutter characteristics or dynamic res- 
ponse of these compressors. Moreover, the basic understanding 
of the flow gained by such a study can lead to better design 
criteria . 

Two basic cascade configurations can be distinguished at 
the outset; i.e., cascades with either subsonic or supersonic 
leading-edge locus. Although the former is the case of pri- 
mary interest in current research, in this paper the case of 
supersonic leading-edge locus is treated in order to form a 
basis for extension to the more complicated case of subsonic 
leading-edge locus. This method of approach also has the 
advantage that the problem of oscillatory wind tunnel inter- 
ference, already analysed in previous work, is contained in 
the theory as a special case. 

The interference problem of linearized supersonic flow 
past airfoils oscillating between solid wind tunnel walls 
was considered by Miles (1956) who derived a solution using 
Laplace transform techniques. Drake (1956) treated this case 
for wind tunnels with free jet boundaries and later Drake 
(1957) gave a solution to this interference problem for 
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porous-wall wind tunnels also using Laplace transform methods. 
In an extension of Miles’ work. Lane (1957) presented a solu- 
tion for supersonic flow past oscillating cascades with super- 
sonic leading-edge locus again using Laplace transform tech- 
niques. Further work in this area was done by Hamamoto (1962) 
Using Teipel’s (1962) method of characteristics approach, 
Platzer and Pierce (1970) made an analysis of oscillatory 
supersonic wind tunnel interference and with the help of the . 
high speed computer were able to predict the pressure distri- 
bution along an airfoil oscillating with arbitrary frequency 
in solid, free jet, or porous-wall wind tunnels. Platzer 
(1971) will present an elementary analysis of porous-wall 
wind tunnel interference effects which generalizes Sauer's 
(1950) solution for a single airfoil oscillating at low fre- 
quency in an unbounded linearized supersonic flow. Platzer 
and Chalkley (1972) further extended this solution to form an 
elementary theory that together with a method of characteris- 
tics procedure (based on Platzer and Pierce, 1970) is used to 
analyse the supersonic flow past oscillating flat plate cas- 
cades having supersonic axial velocities but otherwise arbi- 
trary stagger angle. 

In this thesis the detailed description of both the ele- 
mentary theory and the method of characteristics that was not 
possible in Platzer and Chalkley (1972) is given. With the 
aid of the high speed computer, the method of characteristics 
procedure is used to predict the pressure distributions and 
aerodynamic forces and pitching moments on a typical cascade 
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blade. Both a tvo-degree-of -freedom (torsion and bending) 
flutter analysis and a one-degree-of -freedom (torsion) flut- 
ter analysis are presented, with numerical results given for 



the case of the torsion 
(valid only for low fre 
provide a convenient ch 
as well as to easily pr 
boundaries for a slowly 
interference from super 
and from a larger airfo 
case . 



al flutter. The elementary theory 
quencies of oscillation) is used to 
eck of the characteristics approach 
edict a number of dynamic instability 
oscillating airfoil subjected to 
sonic wind tunnel walls in one case 
il in close proximity in yet another 
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II. PROBLEM FORMULATION 



The cascade is considered to be one of two-dimensional 
flat plates with each blade performing equal low amplitude, 
simple harmonic oscillations of the same mode with the same 
interblade phase angle. As shown in Figure 1, it is aligned 
in the x-y coordinate system such that one blade lies along 
the x-axis with its leading-edge at the origin. 



u 



OO 



y 




\ 

\ 

\ 



Figure 1 



The cascade is further considered to have supersonic 
leading-edge locus and, except for this restriction, arbi- 
trary stagger, such that, 

tan cot a (1) 

With this assumption, the flow between any two adjacent 
blades can be used to describe the flow through the entire 
cascade. Since disturbances cannot travel upstream of the 
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Mach lines from the leading-edge of a blade, interference 
from other than adjacent blades is not possible in the super- 
sonic flow between them. Hence, only the flow between the 
blade along the x-axis and the one above it is considered. 
This flow is assumed to be the non-viscous flow of a perfect 
gas governed by the continuity equation. 



Dp , ,3u . 3v> „ 

D? + p + 37> ■ 0 



( 2 ) 



the Euler equations, 



2P. + I iP. = n 

Dt T p 3x 

+ I Ip. = o 

Dt p 8y 



and the energy equation, 



M = o 

Dt 



(3a) 

(3b) 

(A) 



Because of the assumption of small amplitude oscillations 
of the cascade blades, all flow quantities are considered to 



be small perturbations linearly superimposed on freestream 

quantities. The velocities are written as 

u = u + u (5a) 

00 

c - c + c (5b) 

while the pressure and density perturbations are written as 

Ap - p - (6a) 

Ap = p - Poq (6b) 



With the assumption of isentropic flow, the local velo- 
city of sound can be defined as, 

c 2 = 4^ ' 



dp 



(7) 



12 



further 



= constant 



Taking the total differential. 



( 8 ) 



^7 d P - YP -7 +V dp = 0 

p« p' 



(9) 



or 



hence , 



it = 

dp 



Y 



£ 

P 




( 10 ) 



(ID 



The flow boundary condition along the blade surfaces [ 
shown in Ch. 5, Bisplinghof f , Ashley, Halfman (1955) for a 
simple airfoil] is as follows: If the surface over which 

fluid flows has the equation, 

F (x, y, t) = 0 (12) 

particles in contact with the surface must have the same n 
mal velocity as the surface. Stated differently, the rate 
change of F is zero when the motion of a particular fluid 
element is followed along the surface, or 



DF 

Dt 



= 0 



(13) 



If an arbitrary, thin, two-dimensional cascade blade i 
considered, the equation of the upper surface can be writt 
F u (x,y,t) = y - y u (x,t) = 0 (14) 

where y = 0 is the mean camber line and y is the distance 
from the mean camber line to the upper surface. Likewise, 
the equation of the lower surface can be written. 



as 



the 



or - 
of 
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en. 
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0 



(15) 



F l (x,y,t) = y - y L (x,t) = 



At y = y 



u 



and a t y = y 



DF 

i 

Dt 

L 

Dt 



3y 



u 



3 1 



3 1 



- u 



3y 

J u 



- u 



3x 



3x 



+ v = 0 



(16) 



+ v = 0 



(17) 



since 



3y 



u 






= 1 



3y 3y 

Thus, the normal flow velocity can be written, 

9y, 



(18) 



9y, 

3 t 
3 yi 



U - ' U 

V = Tt~ + U 3x“ 



L - 8y L 
V = — + u 3T- 



at y = y u (x,t) 



at y = y L (x,t) 



(19a) 



(19b) 



By applying the assumption of linear perturbations on the 
freestream, the normal velocity equations become 



3y 



v = 



v = 



u 



3y 



3 1 
3 1 



+ u 



u 



3y 



00 3 x 

3 y T 



+ u 



00 3 x 



+ u 



+ u 



u 



3x 

3x 



at y = y (x,t) 



u 



at y = y L (x,t) 



(20a) 



(20b) 



Since a thin blade is being considered y and y T are small, 

u L 

thus the equations become, cancelling the higher order terms. 



3y 



v = 



V * 



u 



3 t 

9 y 7 



+ u 



3y 



u 



+ u 



w 3x 

9 y 7 



at y = y u (x,t) 



at y = y L (x,t) 



(21a) 

(21b) 



3 1 ’ “co 3t 

This normal velocity can be further written as a Taylor 
series expansion of the normal flow velocity at y = 0. Thus, 
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(22a) 



v(x,y u ,t) = v(x,0 + ,t) t y u — 



+ 'Ll a 2 v(x,o' ,t) 

21 - - T 



and 



ay 



v(x,y L ,t) = v(x,0 _ ,t) + y L 9v |y X ° * ^ 



y 2 9 2 v (x , 0 ,t) 

+ TT 9y 2 + 



(22b) 

Again with the assumption of a thin blade and small linear 



perturbations on the freestream: y u , y^ , and v are small and 
higher order terms may again be cancelled from the equations 
leaving : 



v(x,y u ,t) = v(x,0 + ,t) 




(23a) 


v(x,y L ,t) = v (x , 0~ , t ) 




(23b) 


Thus , 






9y 9y 

v (x , y , t ) - 3t + Uoo 9x 


at y = o"*" 


(24a) 


3y L 3y L 

v(x,y,t) = 3t + Uoo 3x 


at y = 0 


(24b) 


With the assumption of simple harmonic 


motion of 


the blade. 


y u (x,t) = y u ( x ) eiwt 




(25a) 



y L (x,t) = y L (x) e 



ioot 



and the normal velocity equations can then be written, 
v(x,y,t) = [iwy u (x) + u ^ " + 



v (x, y , t ) = [iuy L (x) + u 



9 y 

u ] 
9x J 


iiot 


at y 


e 


3y L 

b i 

co 9x J 


iajt 


at y 


e 



(25b) 

(26a) 

(26b) 



Since the flow is identical between any two adjacent 
blades , 
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at y 



d 



(27) 



v (x , y , t) = 
In plunge. 



9y ] 

Ii03y L (x) + 



J e- 



the lower blade oscillation about y 

, iwt 
y = - h e 
4 o 



0 is , 
(28) 



hence , 

v(x,y,t) = [-itoh Q ] e i ' UJt at y = 0 + 

while the upper blade oscillation about y = d is, 

, i (tot+<$ ) 

y = - h e 
o 

where <5 is the interblade phase angle. Hence, 
v(x,y,t) = [to sin i h Q - iii) cos 6 h Q ] e* wt 



(29) 



(30) 



(31) 



at y = d . 

In pitch, the lower blade oscillation about y = 0 is, 

y = - 9 (x-x ) e iU3t (32) 

o o 

hence , 

v (x , y , t) = [-u 9 - iw0 (x - X )] e iwt (33) 

w o o o 

at y = 0 + , while the upper blade oscillation about y = d is, 

y = - 0 (x-B-x) e i ( wt+ ^) (34) 

J o o 



hence , 

v(x,y,t) = [-0 {u cos 6 -to sin 5 (x-B-x ) } 

- i0 (u sin 6 + w cos S (x-B-x )}] e^ U)t (35) 
o 00 o 

at y = d . 

The pressure distribution on the blades can be written in 
terms of the local sonic velocity as follows: In terms of 

perturbation quantities, Eq. (11) can be written, 
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( 36 ) 



1 + A P 

c 2 + 2 c C = Y — T — 

« P. a + M) 



cancelling the higher order c 2 term. Since ^ is less than 

p 

' OO 

one, Eq. (36) can be expanded in the geometric series, 

1 



= i _ M + Ao =1 _m 

(1 + M) P„ P=o Poo 

pj 



(37) 



neglecting the higher order terms. Substituting this in 
Eq. (36) gives. 



2C C co s T 



(1 - 



A£ 

Y Ap 



) Ap 



(38) 



Since Ap and Ap are small, 

A£ * _i _1_ 

Ap djD -2 

dp 



(39) 



Substituting this in Eq. (38) gives after some algebraic 
manipulation , 



c c (- 

or% ' 



y-1 '*<» v , , 2c 



1 + 



1 N 1 A 

_) = — Ap 

“ no 



(40) 



(Y-Dc 



Since 2c/(y-l)c oo is small, 

the geometric series, 

1 , 2c 



1 + 



2c 



may be expanded in 



(Y-1) c 



1 + 



2c 



= 1 - 



+ Ittt-tt ] 2 - 1 - 2c 



(y~l ) c 



l (Y-l)c 



(Y-Dc 



(41) 



neglecting higher order terms. Hence, Eq. (40) can be writ- 



ten, neglecting higher order terms, 

2 . 

P P co y-1 P°° C °° C oo 



(42) 
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Finally, any non-dimensional quantities so stated in the 
paper are made such with reference to the blade chord length 
and the freestream velocity unless otherwise defined. 
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III. METHOD OF CHARACTERISTICS 



In the method of characteristics it is desirable to as- 
certain the coordinate system across which all possible dis- 
continuities in flow properties may occur. Along this 
coordinate system the equations of motion of the flow field 
can then be treated as ordinary differential equations that 
are solvable by classical or numerical techniques (e.g., 
finite differences). To obtain the equations of this coor- 
dinate system (the characteristic directions) in the (x,y) 
plane, the equations of motion are written in terms of the 
arbitrary intersecting coordinates, 

C = i (x,y) (43a) 

and 

n = n (x,y) (43b) 

If the first derivatives of the dependent variables, u, v, 
and p (or c) with respect to £ are made indeterminate across 
lines of n = constant, and the first derivatives of the de- 
pendent variables with respect to X] are made indeterminate 
across lines of £ = constant, then any possible discontinui- 
ties in the first derivatives will occur across these lines. 
These lines are then the characteristics and their equations 
are obtained in the evaluation of the inde terminacies . 

Consider, .first, two-dimensional, steady flow* The gov- 
erning equations of motion are the continuity equation. 



9 (pu) + Kpv) _ o 
3x 9y 
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the Euler equations, 

3 lH + v 1« + I l£ . o 

3x 3y p 3x 



- 3 v . 3 v , 1 3p 

U + V + — -r*- = 

3x 3y p 3y 

and the energy equation, 

- 3s 



(.45a) 

(45b) 



, 3s f, 

u + v = 0 
3x 3y 

Along a streamline this equation may be written, 



(46) 



~ 3 p , 3 p _, r ~ 3 p , 3 p 

U TT*“ + V -tT- = r 2 [u -r 1 - + v 



} 



(47) 



3x v 3y " c L “ 3x ' v 3y ' 

Substituting in Eq. (44), the continuity equation becomes, 
3u . 3v- 



p tS + - - r 2 If + ^ 



(48) 



In terms of the arbitrary coordinates, Eqs. (43), the 
continuity equation becomes, 



p£ u f + p? v f + — [u£ + vC ] P 

x t. y % -2 x y 



= -pri u - pri v [ u n + v r) ] p 

x n K 'y n -2 x y J F n 

In like manner, the Euler equations become, 



(49) 



p[u£ + v£ ] u ^ + £ p _ = -p [un + vri ] u - p n (50a) 

^x s y £ ?x r ^ x y n 0 x 

p[u£ + v£ ] + £ p^ = -p[uri + vri ] v - p n (50b) 

<^ x s y J r s y £ x y n 0 y 

Eq . (49) and Eqs. (50) form a system of three simultaneous 

equations in u , v , and p . Solving for p by Cramer's 

rule gives a ratio of two determinants. The denominator is 

the determinant of the matrix of coefficients of u^, v^, and 

p^ formed by Eq . (49) and Eqs. (50), while the numerator is 
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this determinant with the right hand side of the equations 
substituted for the coefficients of p . 

5 

Since across lines of r\ - constant, p is indeterminate, 
both the numerator and the denominator must equal zero. Set- 
ting the denominator equal to zero and expanding the deter- 
minant gives. 



+ l^ 2 - c 2 ) £ 2 + 2uv +(v 2 - c 2 )^p = o 



iz = X = t 

dx u 



an 



S = tan ^ 1 a) 







(51) 


quati 


ons of 


all 


y) pi 


ane . 


These are. 






(52a) 






(52b) 


with 


the x 


-axis and 


is tic 


direc 


tions are 


e arb 


i trary 


and in- 


were 


wri 1 1 


en such 


ie unknowns , 


making 


tan t 


would 


result in 


i. Eq. (52 


a) is 


: Eq s . 


(52b) 


are the 


[ach 1 


ines . 


Tradi- 


l = 


constant and 



T1 = constant. 

The compatibility relation which relates changes in u, v, 
and p along the (£,ri) characteristics can be obtained by 
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setting the numerator of p, equal to zero. Expanding the 
determinant and solving the resultant equation yields after 
some algebraic manipulation, 



9u 9 v 7 1 . 9p _ 

v 3W' u 3rr + p co,: “i^ = 0 



(53) 



This derivation is shown in detail in Chapter 8 of Oswatitsch 
(1956) • These relations must be satisfied along the (£,r|) 
characteristics . 

The physical significance of these characteristics is 
that they act as boundaries across which any discontinuities 
will occur. In this flow field where streamlines and Mach 
lines are the characteristics, a discontinuity is readily 
apparent across the Mach line emanating from the leading-edge 
of a slender airfoil or across a slip plane (streamline) 
emanating from a Mach stem. Characteristics can further be 
described as information carriers in that along them all dis- 
turbances propagate and all flow quantities can be determined 
given their values in the domain of dependence. It should 
also be noted that the characteristic directions are com- 
pletely independent of the type of coordinate system used to 
describe the flow. They are determined solely by the physical 
nature of the flow field as expressed in the equations of 
motion . 

The problem of two-dimensional, unsteady flow over a flat 
plate as shown by Teipel (1962) can be considered in the same 
manner as before. The governing equations of motion [Eq. (2), 
Eq. (3), and Eq. (4)] are written in an arbitrary coordinate 
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system across which variation of flow properties are indeter- 
minate. The evaluation of these indeterminacies will again 
give the characteristic directions of these coordinates. 



The • arbitrary coordinates are given in Eqs. (48). Eq.(2), 
Eq. (3), and Eq. (4) can be rewritten in a form more suitable 
for transformation by using the local sonic velocity as an 
unknown. Assuming that air is a perfect gas, Eq. (11) can 
be written, 

(54) 



c 2 = Y — T 

' m. 



where is the universal gas constant and 7U is the molecular 

weight. Taking the total differential and dividing by c 2 
gives , 

2 — = Y — (55) 



d c 


ft 


dT 


- 


Y w 


~ 2 


c 




C 


(54) 


into 


this 


dc 


dT 





(56) 



In like manner, the total differential of Eq. (11) is, 

2c dc = y - Y 

P P P P 



(57) 



and from Eq. (11) again. 



o dc _ dp dp 

- ' p " p 



(58) 



From the First Law of Thermodynamics and the definition of 
entropy, 

(59) 



Tds s du - ^ dp 
P 2 
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Since du = c^ • dt, this may be written as, 



Tds 



c dX - RX ^ 

v p 



(60) 

( 61 ) 



where R is the gas constant (.c -c ) and c and c are the 

p V p v 

specific heats of the gas at constant pressure and volume, 
respectively. Dividing Eq. (61) by Ic^, substituting from 
Eq . (56) and Eq. (59), and then taking the substantial deri- 



vative of 



-S with respect to time gives. 



I ££ = 2y i - (v-i) I = o 

c Dt Y - Dt ^ Y ' p Dt 



(62) 



since entropy is held constant. This can be written substi- 
tuting from Eq. (11), 



Dp = _2_ 



pc 



Dc 

Dt 



(63) 



Dt y-1 

Introducing the small perturbation assumption and substitut- 
ing from Eq . (7) and Eq . (63), Eq. (2) can be written, can- 

celling higher order terms, 



2 3c 2 3 c , 3 u 3v 

y-1 3t y-1 u °° 3x c °° 3x c °° 3y 



= 0 



(64) 



This is the form of the continuity equation that will be used 
in the coordinate transformation. Again dividing Eq. (61) by 
Tc^, substituting from Eq. (56) and Eq . (58), and then taking 

the partial derivative with respect to x (holding entropy 
constant) gives, 

- c ~ C65) 

p 3x y-1 dx 

Substituting in the first Euler equation, Eq. (3) becomes, 
cancelling higher order terms, 



24 



3_u 3_u 

3t u °° 3x y-1 '' 00 3x 



2 c 0 



C66) 



In like manner, the second Euler equation, Eq. (A), becomes, 
3v 



, 3v , 2 3c _ _ 

3t U °° 3x y-1 Cco 3y 



(67) 



The continuity equation and the Euler equations thus form a 
system of three simultaneous equations in u, v, and c. With 
assumption of simple harmonic motion, Teipel (1962) intro- 
duced his amplitude functions. 



U(x,y) e 



rut 



u - u 



u 



( 68 ) 



V (x,y) e 



C (x,y) e 



io)t 



itot 



/M 2 -l u °° 



2 1 c - c 



^2 C oo 



(69) 

(70) 



where U, V, and C are complex. 

Substituting in Eq. (64), Eq. (66), and Eq. (67) gives. 



+ /M 2 -l -r— + M 2 -r— + ikM 2 C = 0 

3x 3y 3x 

3U , 3C .... . 

7T— + x— + ikU = 0 

3x 3x 



3 V 



3 C 



+ ikV = 0 



(71) 

(72) 

(73) 



3x /m 2 -1 3y 

where k is the reduced frequency. 

Transforming this system into the new coordinate system, 
£(x,y) and n(x,y), and considering now a system of three 
simultaneous equations in unknowns, U^, V^, and C^, the equa- 
tions may be written, (74) 

S x U 5 + cot a S y V ? + M 2 C x C 5 - -U n - cot a V^-M 2 ^ " ikM2c 



K u_ 

x C 



+5 c c = -u n -c n - iku 
x c n x n x 



(75) 
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07 * 0 , 



£ x V £ +tan a c y c £ * -V n - ttau a n y - a«K 

As before , and are imade firde term ora £ 3 ^ aa^riKSSS 

lines of r\ * constant. This forces any discontianjJLtiess inr 
U, V, and C to lie across lines af r| — constant. Oiree cnirdi-- 
tion of indeterminacy is that the determinant of thre^ matrix: 
of coefficients of , and in the above ejquattonss be' 

zero. If this determinant is expanded and then set: eaqoxall to: 
zero, the resulting equation is, 

- Sx - e* - k2 ^ ] - 0 

As before the solutions to this equation give the thneef ckar:- 
acteristic directions in the physical (x,y) plane: 

= 0 (738a} 

4^- = + tan. a (738 b.) 

ox — 

These coincide with the results af the steady flow. anaiy - 
sis when it is remembered that the small perturbation assump- 
tion forces the streamlines to be paraJLlel with tha x— axis 

U = 0 ). 

In order to find the compatibility relations of the first 
derivatives of U, V, and C along the characteristics, the 
equations of motion are written along lines of C = constant, 

T) = constant, and along streamlines, that is, for 



str(x,y) = constant. 


lx = 

3-x 


0 


(79 a) 


£ (x,y) = constant. 


lx - 


1 


(79h.) 


3-x 


»<ll z -l 
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r\(x,y) = constant. 



(79c) 



Iz « _ 

9x 



/M 2 -l 



If an arbitrary function of x and y in the equations is 



3 f 

denoted by f then, holding str constant becomes, 



(M) 

^ 9x^ s tr 



H 

9x 



(80) 



9f 



Likewise, holding £ constant becomes 

fM') = 11 . 1 M 

9x 5 9x * 9y 

3 f 

and, holding n constant becomes. 



(81) 



(11) 9 f 



9 f 



9x n 9x 9y 

Combining these equations and rewriting, 

9 f = ,9f ' 

9x *9x str 



11 = lr (II) + ( 11 ) ] 

9x 2 l ''9x'£ '■9x ; iv l 



(82) 



9 f 



. 9 f . 



. 9 f . 



fy - i ^ 2 - 1 ~ 



(83a) 

(83b) 

(83c) 



Substituting these relations in the equations of motion gives. 



tk!?>, + (S.i + t(m 2 - imS, - (S.i 



9x £ 



+ 7 M 2 [(|§)^ + (|^) n ] + ikM 2 C = 0 



9x T) 2 

2 M * [ v 9x^ 

1 

TV 



' 9x £ 



9x T) 



(84a) 



I + ^n 1 + §«<?£>? + + iku - 0 < 84b > 



+ + + ikv ‘ 0 < 84 ‘> 

Subtracting Eq. (84b) from Eq. (84a) gives an equation that 
is first added to Eq. (84c) and then subtracted from Eq . (84c). 

’ The result is two compatibility relations that must be 
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satisfied along lines of constant £ and lines of constant r\ . 
These are, 

(M 2 C - U)] = 0 



(—) + C— ) + ik rv + — - — 

( 3x^ + C 3x } £ + ik IV + 



a 8V\ /• 3 C ^ rTT 1 

' c 9^n + lk Iv .' ~ 



(M 2 C - U)J - 0 



(85a) 

(85b) 



M -1 

The third compatibility relation is obtained by considering 
Eq. (84b) along a streamline, 

O n 

( 86 ) 



Ostx + <!£>.« + ikD - 0 



If irrotationality is assumed, the equation of irrotationali- 

ty. 

3u 3v 



= 0 



(87) 



3y 3x 

may be substituted for the second Euler equation in the sys- 



3U 

3y 



as 



H + /FT! |I + m 2 

3x 3y 



on of 


simple h 


armonic 


motion, i 


ipel 1 s 


ampli tu 


de func 


t ions , 


- i 


F - 0 

dx 




(88) 


equat 


ions can 


then b 


e written 


i£ + 

3x 


ikM 2 C = 


0 


(89a) 


■ + i£ 

3x 


+ ikU = 


0 


(89b) 


/M 2 - 


■ i H . 

3x 


0 


(89c) 


the arbitrary 


coordinate system 



3U 

3y 



£ an ^ n, as before, the system can be considered one of simul- 
taneous equations of U , V , and C .. By making any of these 

£ 5 £ 
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derivatives indeterminate the results of setting the deter- 
minant of the matrix of their coefficients equal to zero are 
Eqs. (78), the characteristic directions. 

By writing the equations of motion, as before, along 
lines of £ - constant and T| - constant and performing the 
same algebraic manipulations, the two compatibility relations 
can be determined: 



- 




+ ik 


' M 2 


(U - 


C) = 0 


(90a) 


T— ( 

I 

CM 

a 






+ ik 


M 2 


(U - 


c)=o 


(90b) 


H 2 -l 



The third compatibility relation is, as before. 



(P-) + (|£) . + ikU = 0 

3x str dx str 

The equations of motion are now reduced to 
taining only derivatives with respect to x. It 
possible to solve these equations using finite 
in the following manner: Separating the real f 

nary components, the equations can be written. 



( 86 ) 

a system con- 
is, thus, 
differences 
rom the imagi- 



aU R dC Tt 

(- 3 -*) „ + ( 3 — O „ - kU T 

dxstr 3 xstr I 



= 0 



3U 9 C 

^str + ^.tr + kU R 



= 0 






au R 3v r 

( 7 — O - (- 5 -^) - k 

9x K 9x 5 M 2 -l 



M' 



3u av T 

C a~) - + k 

9x 5 9x 5 M 2 -1 



aU R 3V R 

( 3^n + <T5T>n " k .a 






M -1 



<V V 



<V V 



cu x- V 



= 0 



- 0 



= 0 



(91a) 

(91b) 

(91c) 

(91d) 

(91e) 
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(91f) 



au T av T 
+ C 5T ) - 



+ k 



M 

n z -i 



<V V = 



To write these equations in finite difference form, a 
computational molecule as shown in Figure 2 will be used. 
With the velocity values at ^12* anc ^ ^21 known, a H tlle 

values at P ^2 can be calculated. The distance Ax is arbi- 
trary . 




Figure 2 



If F is an arbitrary flow quantity, its partial deriva- 



tives can thus be written, 

(— ) 

s tr 

ell) 



(H) 

'•ax'n 



F - F 
22 *11 

Ax 

F - F 
22 *12 

yAx 

F 22 " F 21 



(92a) 

(92b) 

(92c) 



The values of the flow quantities in the equations are aver- 
ages, such that, 

(F) str = 2 ^ F ll + F 2 2 ^ (93a) 
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(93b ) 



00 
■ 5 



00 



n 



2 ^12 + F 22 ^ 



2 ^21 + F 22 ) 



(93c) 



With these relations substituted into the equations, the sys- 
tem becomes, after some algebraic manipulation. 



where : 



and , 



22R 


+ 


C 22R 


— 


A I U 22I 


= 


K 12R 




(94a) 


221 


+ 


C 22I 


+ 


A I U 22R 


= 


K 12I 




(94b) 


22R 


- 


V 

V 2 2R 


- 


B I U 22I 


+ 


B I C 22I = 


K 34R 


(94c) 


221 


- 


V 

221 


+ 


B I U 22R 


- 


B I C 22R = 


K 34I 


(94d) 


22R 


+ 


v 

22R 


- 


B I U 22I 


+ 


B I C 22I = 


K 56R 


(94e) 


221 


+ 


v 

221 


+ 


B I U 22R 


- 


B I C 22R = 


K 56I 


(94f) 



12R 


— 


U 11R 


+ 


C 11R 


+ 


A I U 11I 




(95a) 


121 


= 


U 11I 


+ 


C 11I 


- 


A I D 11R 




(95b) 


34R 


= 


U 12R 




V 12R 


+ 


B I (U 12I- 


C 12I ) 


(95c) 


341 


= 


U 12I 


- 


Vl21 


- 


B X (D 12R- 


C 12R ) 


(95d) 


56R 


= 


U 21R 


+ 


V 21R 


+ 


B I (U 21I- 


C 2H> 


(95e) 


561 


= 


U 21I 


+ 


V 21I 


- 


V U 21R- 


C 21R^ 


(95f) 



Aj = y k Ax 

- X * — 

1 4 M*-l 



Ax 



(96a) 

(96b) 



As shown in Teipel (1962), solving these equations gives, 
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u 



a-A I B I )[i(K 34R+K56R )-B lKl 2 I 3+2B I [|cK3 4I+ K 5 6 I )+B I K 12R ] 



22R 



+ C2B I )2 



(97 a) 



U 



^ 1 “ A I B I ) *2 ( - K 34l +K 5 6l )+B I K 12R'* _2B I^2 ^ K 34R +K 5 6R^ " B I K 12I^ 



221 



(1-AjBj.) 2 + (2B r ) 2 



V 22R " 2 (K 56r" K 34R> 



V 22I 2 ^ K 56I~ K 34I^ 



(97b) 
(97c) 
(9 7d) 



C 22R K 12R U 22R + A I U 22I 
C 22I = K 12l“ U 22l“ A I U 22R 



( 97e ) 
( 9 7 f ) 



These are the finite difference equations for a point in the 
general flow field. 

At a point on the top of a cascade blade or solid bound- 
ary, the normal velocity V ^ is prescribed by the movement 
of the blade. This velocity was given in the Problem Formu- 
lation. Here the computational molecule is shown in Figure 3. 




Figure 3 
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The applicable equations are, 



22R 


+ 


C 22R 


A I U 22I 


K 12R 




(98a) 


221 


+ 


C 22I + 


A I U 22R 


K 12I 




(98b) 


22R 


- 


B I U 22I 


+ B I C 22I 


= K 56R- 


K 34R 


(98c) 


221 


+ 


B I U 22R 


B I C 22R 


= K 56I- 


K 34I 


(98d) 



where an d K,-^ are as before, but. 



K = V 
34R 22R 



K = V 
341 221 



(99a) 

(99b) 



as given by the boundary conditions. Solving for these 
equations gives, 

(!- A I B l) [ K 56 r _K 34r" B i K 12I ■' + 2B I ^ K 56l _K 34l +B I K 12R-* 



U 



22R 



U 



221 ' 



(1-A I B I ) 2 + (2B i ) 2 

( 1 ~ A I B I ) t K 56l~ K 34l +B I K 12R-* ~ 2B I ^ K 56R~ K 34R~ B 1 K 12I-* 
(1-A I B I ) 2 + (2B ] .) 2 



C 22R K 12R U 22R + A I U 22I 



(100a) 

(100b) 

(100c) 



C 22I K 12I U 22I A I U 22R 



(lOOd) 



^22R an< * ^221 are known f rom the boundary conditions. These 
are the finite difference equations for the flow velocities 
at the top of a solid boundary or cascade blade. 

Likewise, at the bottom of a solid boundary or cascade 
blade, the normal velocity, V i s prescribed by the movement 
of the blade, and was given in the Problem Formulation. The 
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computational molecule, different from that at the top of a 
blade, is shown in Figure 4. 




Figure 4 



Here, the applicable equations are, 



U 22R 


+ 


C 22R “ A I U 22I 


K 12R 


(101a) 


U 2 21 


+ 


C 22I “ A I U 22R 


K 12I 


(101b) 


U 22R 


- 


B I U 22I + B I C 22I = 


K 56R + K 34R 


(101c) 


U 22I 


+ 


B I U 22R _ B I C 22R = 


K 56I + K 34I 


(lOld) 



where and are as originally stated, but, 

K 56R - V 22R < 102a > 

K 56I ‘ V 22I < 102b > 

as given by the boundary conditions. Thus, the finite differ- 
ence equations for the flow velocities at the bottom of a 



solid boundary or cascade blade are. 



U 



22r' 



( 1 - A I B I )[K 56R H-K34 R -BiKi2i3 + 2 Bi[K 5 6i + K34i + B iKi2 R ] 

Cl-A ] .B I ) 2 + (2B X ) 2 (103a) 
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(103b) 



U 



a - A X B I ) I K 56I +K 34l +B X K 12K J - 2B X lK 56 R +K 34R' B I K 12^ 



22X 



(1-A^) 2 + (2B X ) 2 



C 22R" K X2R " U 22R + A I U 22I 



C 22I K X2X " U 22X ” A X U 22R 



Cl03c) 

(103d) 



V 00 and are given from the boundary condition, 

Z Z K Z Z 1 



To obtain the conditions on the initial left-running Mach 
line, assume that and as shown in Figure 2 are just 



n = K cz: = 0 and A_ and B T go to zero as Ax goes to zero, 

i Z 3 O 11 



in the freestream and then let Ax shrink to zero. Since 
K ] 

the initial finite difference equations for £ = constant are, 

(104a) 



U 2 2 + C 22 * 0 



Thus , 



U 22 ~ V 22 + V 12 " ° 



U 22 + V 22 " ° 



U 22 V 22 ” C 22 



(104b) 

(104c) 

(105) 



Eq. (90a) can then be written. 



,3U. 

- lku 



Integrating gives. 



Thus , 



U = U 



exp [ -ik 
x=0 



M" 



M 2 -1 



x] 



(106) 



(107) 



U 22R " ~ V 22R^°^ COS % 2 



M' 



M - 1 



x)- V 22I C0)sin (k — 



M' 



x) (108a) 



H - 1 



U 22I = -V 22l (0)cos(k 






M - 1 



x)+ V 0 op-CO) sin (k 

22K. 2 






x) (108b) 



r- l 
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(108c) 



V = - u 

V 22R 22R 



V 22X ~ " U 22I 



C 22R " " U 22R 



r = _ tt 

^221 U 22I 



(108d) 
(108e) 
(108f ) 



For the initial right-running Mach line, and 

Figure 2 are just in the freestream and Ax is brought to zero. 



Since K, „ = K.. = 0, and A T = B_ = 0, 

12 34 II 



°22 + C 22 * 0 



Thus , 



U 22 - V 22 - 0 



U 2 2 - U 21 + V 22 = ° 



U 2 2 V 22 “ °22 



Eq. (90b) can then be written. 






Integrating gives. 
Thus , 

U 22R = V 22R (0)cos(k 



U = V(0) e t_lk „2 



M- 



x] 



M z - 1 



M' 



M 2 - 1 



U 22I = V 22I (0)cos(k — 



M' 



M - 1 



(109a) 

(109b) 

(109c) 

( 110 ) 



( 111 ) 



( 112 ) 



x)+V„ (0)sm(k 

Z Z 1 _ ,2 


x) 


(113a) 


M - 1 






x) " V 22R C ° )sin(k ~~~ 


x) 


(113b) 


M 2 - 1 






22R = U 22R 




(113c) 


22R =_U 22R 




(113d) 
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(113e ) 



V 22 X " U 22I 
C 22I " " U 22I (113f ) 



The pressure distribution along the surfaces of the blades 
can be determined from Eq. (42). Dividing both sides by 
and using Teipel's amplitude function for C, the equation 
can be written, 



P (x,y) 



yM 2 C (x,y) 



x=0 

x=0 



(114) 



where 



P (x,y) e 



icot 



P " P r 



Thus the 
top surface 
be computed 
along lines 



non-dimensional pressure distribution along 
and the bottom surface of the cascade blade 
in a point by point calculation of U, V, and 
of £ = constant and f| - constant. 



the 

can 

C 
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IV. ELEMENTARY THEORY FOR SLOWLY OSCILLATING CASCADES 



Sauer's (1950) solution for an airfoil 
low frequencies in an unbounded supersonic 
to the flat plate cascade to form a theory 
parison with the method of characteristics, 
rick (1957), the equations of motion of sup 
a flat plate oscillating at small amplitude 
assuming irrotational flow, as 

cot 2 a$ - $ + 2 - - $ =0 

xx yy c^ xt c 2 tt 



oscillating at 
flow can be applied 
suitable for com- 
As shown in Gar- 
ersonic flow over 
can be written, 

(115) 



With the assumption of simple harmonic motion, 

$(x,y,t) = 4>(x,y,k)e ikt 

and the potential equation becomes. 



(116) 



cot 2 a<b - d> + 2ikM 2 d> - k 2 M 2 <t> = 0 
xx yy x 

where $ is complex. 

If the frequency of oscillation is 

perturbation potential amplitude may be 

series, such that 

<Kx,y,k) = x( x »y) + k 

neglecting higher order terms. Eq. (11 

two simultaneous equations since k cann 

tion, it being independent of k. These 

cot 2 a X "X = 0 
xx ^yy 

cot 2 a ip - X p = -2iM 2 X 

r xx r yy A x 



(117) 



sufficiently low, the 
expanded in a Taylor 

(x,y) (118) 

7) then splits into 
ot appear in any rela- 
equations are: 

(119a) 

(119b) 
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Sauer (1950) showed that a general solution to these two 
equations is, 

X 55 sCz) 



« h(Z) - i 



M 



cos a 



yg(z) 



( 120 ) 
( 121 ) 

( 122 ) 

(123) 

(124) 

where g and h are arbitrary functions of position and are 
equal to zero for z ^ 0 (x,y, z and d are non-dimensional). 

To apply this general solution to the cascade with super- 



z 


s= 


X 


- y 


cot 


a 


h(Z) 


+ 


i 


M 




(z) 


COS 


a yg 


Z 


= 


X 


+ y 


cot 


a 



sonic 


leading-e 


dge locus 


> 


the 


flow f 


ield is divided 


into 


separate 


zones 


as shown 


in 


Fig 


ure 5 . 


The number of 


these 


zones 


is 


dependent on th 


e 


size 


of A . 


For A > 1, there is 


one zone 


along 


each blad 


e; 


for 


0.5 < 


A ^ 1, there are two 


for 0 


.33 


< A 


0.5, ther 


e 


are 


three ; 


for 0.25 < A ^ 


0.33, 


there 


are 


f our ; 


etc. 













If oscillation in pitch only is analysed, the boundary 
conditions on the blade surfaces can be written, in non- 
dimensional form, from Eq. (33) and Eq . (35), 



4> y = -1 - ik (* “ X D ) 

<P = e i<S [-l -ik(x - x )] 



at y = 0 
at y = d 



(125) 

(126) 



where 0 is one. 
0 
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Figure 5 



In Zone I there is no interference from the upper blade, 
hence, Sauer's solution for the single airfoil applies. The 
boundary condition, however, must be satisfied at y = 0, 
hence , 

Xy - - cot a g 0 '(Z) = -1 (127a) 

= - cot a h 1 ( Z) -i — - — g ( Z) =-i (x-x ) (127b) 
r y o cos a & o o 

This gives, 

g Q (Z) = Z tan a (128a) 

h (Z) =-i tan aZ (x + \ tan 2 a) (128b) 
o o l 

where , 

* = g o (z)+ k ^ 0 CZ)-i y g 0 (Z)} ( 129 ) 

In terms of Teipel's amplitude functions the flow properties 
can be related to those obtained in Method of Characteristics 

by, 
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